Nearly frozen Coulomb liquids 
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We show that very long-range repulsive interactions of a generalized Coulomb- like form V[R) ~ 
R~°' , with a < d ( d- dimensionality ) , typically introduce very strong frustration, resulting in extreme 
fragility of the charge-ordered state. An "almost frozen" liquid then survives in a broad dynamical 
range above the (very low) melting temperature Tc which is proportional to a. This "pseudogap" 
phase is characterized by unusual insulating-like, but very weakly temperature dependent transport, 
similar to experimental findings in certain low carrier density systems. 

PACS numbers: 71.30.+h,71.27.+a 



I. INTRODUCTION 

In designing novel materials, lightly doping a parent 
insulator is typically the method of choice. An espe- 
cially intriguing situation is found in ultra-clean samples 
at finite doping, where neither the Andersoni (disorder- 
driven) nor the Mott^ (magnetism-driven) route for lo- 
calization can straightforwardly succeed in trapping the 
electrons. The tendency for charge ordering (CO) then 
emerges as the dominant mechanism that limits the 
electronic mobility. As first noted in early works by 
Wigner— and Mott^, this is precisely where the incipi- 
ent breakdown of screening reveals the long-range na- 
ture of the Coulomb interactions. The corresponding CO 
states proved to be of extraordinary fragility, restrict- 
ing the insulating behavior to extremely low densities 
and/or temperatures^. A broad range of parameters then 
emerges where puzzling "bad insulator" transport char- 
acterizes such nearly-frozen Coulomb liquids. 

Unusual "bad-insulator" transport behavior has 
been observed in many systems. Examples range 
from high mobility two-dimensional electron sys- 
tems in semiconductors)^ to lightly-doped cuprates;^ 
manganitesj^ and even to the behavior of lodestone (mag- 
netite) above the Verwey transition)^ In all these cases, 
a broad range of temperatures has been observed, where 
the resistivity rises at low temperatures, but it does so 
with surprisingly weak temperature dependence. In con- 
trast to conventional insulators, where the familiar acti- 
vated transport reflects a gap for charge excitations, the 
"bad insulator" behavior has been interpreted^ as a pre- 
cursor to charge ordering, leading to very gradual open- 
ing of a soft pseudogap in the excitation spectrum. 

The physical picture of a nearly-frozen Coulomb liq- 
uid has been proposed-on a heuristic level-by several 
authors^"— providing a plausible and appealing inter- 
pretation of many experiments. The interplay of spins 
and charge degrees of freedom in pseudogap formation is 
still a controversial and unresolved problem. Therefore, 
to focus on the corresponding role of charge fluctuations, 
we deliberately ignore any spin effects, and consider a 
class of models of spinless electrons interacting through 
long-range interactions. 
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Figure 1: (Color online) Phase diagram of the half-filled clas- 
sical d = 3 lattice model with interactions V{R) — R~" . The 
charge ordering temperature Tc{a) ~ q, as obtained from 
EDMFT theory (full line) and Monte-Carlo simulations (open 
symbols). The pseudogap temperature T* (dashed line) re- 
mains finite as a — 0; a broad pseudogap phase emerges at 
a < d. We also show T^^ ~ 1 for the same model with 
short-range interactions (dotted line), and T^^^ (dot-dashed 
line) from the classical limit of RPA. The inset shows the cor- 
responding plasmon mode spectral density, which assumes a 
scaling form for a <C 1. The fluctuations of these very soft 
"sheer plasmons" lead to the dramatic decrease of Tc- 



We present the simplest consistent theory of this 
strongly coupled liquid state. We demonstrate that the 
existence of such an intermediate liquid regime ,which 
emerges at /cbTc < ksT ^ Ec (see below), is a very gen- 
eral phenomenon reflecting strong frustration produced 
by long-range interactions. It holds for any interaction 
of the form V{R) ^ R^", both in continuum and lat- 
tice models at any dimension d > 2, with a <^ d. Ours 
is a microscopic theory that substantiates this physical 
picture^^iii based on quantitative and controlled model 
calculations. We present a physically transparent analyt- 
ical description using extended dynamical mean-field the- 
ory (EDMFT) to accurately describe the collective charge 
fluctuations, and benchmark our result using Monte- 



Carlo MC simulations. 



II. OUR MODEL AND THE EDMFT 
APPROACH 

It has long been appreciated^ii^ii^ that in Coulomb 
systems, the CO temperature scale Tc is generally very 



small as compared to the Coulomb energy 



(a being typical inter-particle spacing), which we use as 
our energy unit. For example, for classical particles on a 
half- filled hypercubic lattice ~ 0.1,— while in the con- 
tinuum and classical Wigner crystal Tc ~ 0.01—; similar 
results are obtained both in d = 2 and in d = 3. Such 
large values of the "Ramirez index"— / = Ec/Tc sug- 
gest that geometric frustration plays a significant role, 
reflecting the long-range nature of the Coulomb force. 

To clarify this behavior, we control the amount of frus- 
tration by introducing generalized Coulomb interactions 
of the form V{R)/Ec ~ (R/a)^"'. We consider a lattice 
model of spinless electrons given by the Hamiltonian 

H = - Y,t^AcJ + lY.^{R^,){n^-{n))in,-{n)). (1) 



Here c| and a are the electron creation and annihilation 
operators, tij are the hopping matrix elements, = cjc^, 
and Rij is the distance between lattice sites i and j ex- 
pressed in the units of the lattice spacing. The origin of 
frustration is then easily understood by noting that in 
the classical limit our lattice gas model {rii ~ 0, 1) maps 
onto an Ising antiferromagnet (Si = ±1) with long-range 
interactions. Here, the maximum level of frustration is 
achieved for infinite range interactions (a — 0), and any 
finite temperature ordering is completely suppressed. 

A controlled theoretical approach to our problem is 
available for very long-range interactions [a <ti 1), which 
effectively corresponds to a very large coordination num- 
ber. In this limit the spatial correlations assume a sim- 
plified form 
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iuJn - Efc - S(itJ„) ' 



+ Vk • 



(2) 



where the momentum dependence of the (fermionic) self- 
energy S(iw„) and the irreducible polarization operator 
Il(ifln) can be ignored^^. A conserving approximation 
that formally sums all the corresponding Feynman dia- 
grams is given by the so-called EDMFT formulation,—"— 
where the relevant (local) quantities are computed from 
an auxiliary local effective action 
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eff 
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+ - I dTdT'6n{T)n-\T-T')Snir'), (3) 
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Figure 2: (Color online) (a) Density of states p{iu, T) obtained 
with three different methods: EDMFT (full line), MC (dashed 
line) and SCGA (dot-dashed line). Results are shown for the 
d = 3 half-filled cubic lattice with t — 0, a = 0.3, and two 
temperatures: T = 0.03 ^ and T = 0.25 ^ T*. p{lj,T) 
obtained from EDMFT (full line) agree well with MC results 
(dashed line), while SCGA (dot-dashed line) does not account 
for the pseudogap formation, (b) Both EDMFT (full line) 
and MC results (open symbols) p(lj — 0, T) show pseudogap 
opening (dramatic DOS decrease) at T < T* in contrast to 
SCGA results (dot-dashed line). 



where Qq^^iu}) = iw — A(iw) and Snlr) = n(T) — (n). 
The dynamical effective-medium (EM) functions A and 
Hq ^ represent the respective fermionic and bosonic baths 
coupled to the given lattice site. For a given bath, the 
(local) Dyson's equations stipulate that S = Q'^^ — Gj~J 
and H"^ = Hj^J— H^ ^, where Gioc and Hioc are calculated 
directly from 5cff.The self-consistency loop is then closed 
by relating the local and the EM correlators, viz, Gioc = 
Efc Gfc(S) and Hioc = Efc nfe(ri). 



III. CLASSICAL LIMIT 

The most stringent test for the accuracy of EDMFT 
is provided by examining the classical limit (t = 0), 
where pseudogap formation is most pronounced. Here, 
the EDMFT equations can be solved in closed form,— 
since the "memory kernel" HQ"^(r — r') becomes a time- 
independent constant, Hq^ = I?//?^, and the corre- 
sponding mode-coupling term in Eq. (3) can be de- 
coupled by a static Hubbard-Stratonovich transforma- 
tion. The density correlator then assumes the form 
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Hfc = (4 + D- 
reduces to 



-(3Vk) ^, and the self-consistency condition 



dermis) (4 + D 



(4) 



where we introduced the (classical) plasmon-mode spec- 
tral density (^(e) = J2k ^(^ ^ ^k)- The CO critical tem- 
perature Tc{a) is identified by the vanishing of II^T^ at 
the corresponding ordering wave vector k = Q . The 
mechanism for Tc depression is then easily understood by 
noting that for a ^ 1 the spectral density 1^(8) assumes 

the scaling form i^(e) = a~'^i'{{e—eo)/ct), where Eq « —1; 
the explicit form of the scaling function t'(£ — Eq) corre- 
sponding to the half-filled cubic lattice is shown in the 
inset of Figure 1. It features a sharp low-energy spectral 
peak of the usual dispersive form (^(e) ^ si^--^)/^ only at 
(e — So) < £*{ct), i.e below a characteristic energy scale 
e*{a) ~ a and a long high-energy tail of the form z/(e) ~ 
e^^. Physically, these low energy excitations correspond 
to "sheer" plasmon modes with wave vectorfc Q; the 
scale e*{a) ^ a thus plays a role of an effective Debye 
temperature. Its smallness sets the scale for the ordering 
temperature Tc{a) = a J de i^{s)/e ~ e*(Q^)) in agreement 
with an estimate based on a Lindemann criterion applied 
to the sheer mode^l. 

In the classical limit, the single particle density of 
states (DOS) p(a;,r) = — ImG(a; + iO+)/7r assumes a 
simple bimodal form: 
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with the self-consistently determined parameter D{T) 
setting the scale of the Coulomb pseudogap ("plasma 
dip") Egap = D//3, which starts to open at the crossover 
temperature T* = D/Aj3. We stress that, in contrast to 
the ordering temperature Tc ~ a, both Egap and T* re- 
main finite for a <C 1, since D{T) ?s /3 in this limit. This 
leads to the emergence of a broad pseudogap regime for 
a < d, independent of the precise form or the filling of the 
lattice. Remarkably, since D{T) remains finite as a — ?> 0, 
both the density of states p{uj,T) and the conductivity 
cr(T) (see below) display only very weak a-dependence, 
in contrast to Tc{a) ^ a. 

We benchmark these analytical predictions against MC 
simulations which used careful finite-size scaling anal- 
ysis and (generalized) Ewald summation techniques to 
account for long-range interactions (the detail is in the 
appendix). It was found that EDMFT captures all qual- 
itative and even quantitative features of the pseudogap 
regime for several different values of the exponent a, both 
in dimensions d ~ 2 and in d = 3. The detailed compar- 
ison of EDMFT and MC results will be presented else- 
where; here we illustrate these findings for a c? = 3 half- 
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Figure 3: (Color online) (a) EDMFT phase diagram for a 
half-filled cubic lattice with q = 0.3 as a function of tem- 
perature T and the electron's Fermi energy Ef ~ t. Our 
semiclassical solution is valid above the CO freezing temper- 
ature Tc{Ef) (full line), and the Fermi liquid crossover tem- 
perature Tcvos{Ef) (dot-dashed line). At intermediate tem- 
peratures Tc < T < T* we find well developed pseudogap 
behavior, where transport assumes insulating-like but very 
weak temperature dependence [as shown in (b)]; in the CO 
phase {T <Tc) transport assumes the conventional activated 
form (not shown), (b) Temperature dependence of the con- 
ductivity in the semiclassical regime ( Ef <S 1), where only 
the prefactor CTo = f-^Hj- displays significant f-dependence; 
EDMFT results (full line) again show remarkable agreement 
with results obtained by calculating p{e,uj) in the classical 
limit using MC simulations (symbols). 



filled cubic lattice. Figure I shows how EDMFT accu- 
rately captures the a-dependence of Tc, which is found to 
decrease in a roughly linear fashion as a — ^ 0, while the 
T* 0.25 remains finite, producing a large separation 
of energy scales and a well-developed pseudogap regime. 
Note that the familiar Coulomb interaction {a = 1) lies 
well within the small-a regime. This observation makes 
it clear why our EDMFT theory remains very accurate 
(as noted in previous work^^) not only for a <C 1, but 
also for the physically relevant Coulomb case a = 1. 



IV. GAUSSIAN THEORIES DO NOT CAPTURE 
PSEUDOGAP FORMATION 

The excellent comparison between EDMFT and MC 
results for the DOS is shown for a = 0.3 in Figure 2(a). 
In contrast, the conventional approaches^, which typi- 
cally assume Gaussian statistics for the collective charge 
fluctuations, fail to capture the pseudogap opening at 
T > Tc- For example, the familiar self-consistent Gaus- 
sian approximation ("spherical model"), while predicting 
the exact same Tc as EDMFT, produces Gaussian-shaped 
DOS at any T > Tc, in contrast with MC findings; 
these shortcomings are especially dramatic for a ^ d 
(see Figure 2). The popular "random-phase approxima- 
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tion" (RPA)fi^ which amounts to a non-self-consistent 
Gaussian approximation (SCGA), proves even less reh- 
able in this regime. It grossly overestimates the freezing 
temperature Tc, which is found [dashed line in Figure 2 
(b)] to remain finite even as a — > 0, completely missing 
the pseudogap regime (shaded area in Figure 1). Physi- 
cally, the RPA (Stoner-like) freezing criterion reduces to 
the simplistic Hartree (static mean-field) approximation, 
which ignores the dramatic fluctuation effects of the soft 
collective (sheer plasmon) modes. 



V. BAD-INSULATOR TRANSPORT IN THE 
SEMICLASSICAL REGIME 

We expect the "bad insulator" transport to be best 
pronounced in the semiclassical regime i ^ 1, where the 
Coulomb energy represents the largest energy scale in 
the problem. Here, the pseudogap phase is reached by 
thermally melting the CO state at T > Tc{t). While our 
EDMFT equations are difficult to solve in general, in 
this incoherent regime it is well justified to utilize an 
adiabatic ("static") approximatioufi^ which ignores the 
time dependence of the collective mode. The EDMFT 
equations can then be solved in a manner similar to that 
in the strict classical limit (see above), and we find 



where Po(e) is the bare single-electron density of states 
and A{e,uj) = -ilm(w iO+ - e - E(w iO+))"\ In 
this adiabatic approximation, we calculate conductivity 
in the leading order of t^, in terms of quantities for t ~ 
(/9(e) and A(e, uj)). 

These equations are easy to solve for arbitrary parame- 
ters of our model, but we illustrate our findings in Figure 
3; by showing explicit results for half-filled cubic lattice 
with a = 0.3. Our semiclassical solution is found to 
be valid in a broad pseudogap regime Tc < T < T* , 
which spans almost an order of magnitude in tempera- 
tiu-e (for Ep <S: 1 we find Tc « 0.03 and T* sa 0.25). 
Here the conductivity displays unusual, insulating-like 
[da{T)/dT > 0], but rather weak (almost linear) tem- 
perature dependence [shown in Figure 3(b)[, surprisingly 
similar to that observed in magnetite above the Verwey 
transition. Our microscopic theory confirms the heuristic 
picture first proposed in early work of Motti^ 



VI. CONCLUSIONS 



H(irj) = r^ j d<j)G4iuj + in)G4iu})P{(t)), 



P{<P) = |cxp ( ~^H4> + 1) - ^ln(G4*L.)) 
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G7i(iw) = G^\iuj) + (j)TD. 



(6) 



Physically, the electrons travel in the presence of a static, 
but spatially fluctuating random field representing the 
collective mode. Its probability distribution P{4>) as- 
sumes a strongly non-Gaussian character, reflecting the 
charge discreteness captured by EDMFT, but ignored by 
conventional Gaussian theories such as RPA. 

The semiclassical approximation remains validi^ as 
long as the time-dependence of the density correlator 
n(r) can be ignored, corresponding to 



|(n(o)-n(/3/2))/H(o)|«i. 



(7) 



This criterion provides an estimate for the crossover tem- 
perature Tcrosi below which we expect (at large t) a grad- 
ual crossover towards Fermi liquid behavior. The result- 
ing phase diagram is shown on Figure 3 (a). 

To calculate transport, we use the Kubo formula for 
the resistivity, which within the EDMFT theory assumes 
the form^S: 
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(8) 



We argued that pseudogap behavior in Coulomb sys- 
tems directly reflects strong frustration found in any 
system with very long-range repulsive interactions. We 
demonstrated that a quantitatively accurate strong- 
coupling description of this regime is possible using the 
interaction power a as a small parameter in the theory. 
The corresponding EDMFT equations were solved in the 
semiclassical regime where the pseudogap phenomena are 
most pronounced, explaining "bad-insulator" transport 
found in many puzzling experiments. It should be noted 
that, using appropriately formulated quantum impurity 
solverS)^ the same formulation could be extended to in- 
vestigate low-temperature quantum critical behavior for 
the same class of models. This fascinating direction re- 
mains a challenge for future work. 
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VIII. APPENDIX 



A. Ewald Potential 



To carry out a careful finite-size scaling analysis of the 
DOS, we perforin a two-Gaussian fit 



In order to compute the effective potential of long- 
range interaction 1 / 1 | " in hypercubic lattice 



1 



(9) 



we use an Ewald-type summationS^ with the help of the 
integral representation o^i^i 
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f t'^-^e-''''dt+ [ it-ie-'-'*dt. (10) 
a/2) J J 



where r{a/2) is Gamma function. We switch the first 
term of the integral to a momentum sum because the 
sum does not converge rapidly in real space. Next, we 
use the representation^^ 



J d'^r"^S{r- [Lr 



,])/(^1 



iGi.r 



-'5(r)]/(r0, 



fll) 



where /(r) is any arbitrary function and on the right- 
hand side the summation is over the vectors of the re- 
ciprocal lattice. We then integrate r out and change the 
variable of the integration in the first term t 1/t. The 
final expression of the potential takes the form 



+ 



-ik.f 



k^O 



2e" 
f(f)' 



(12) 



where in each component vector ki = 2-Kni and nitL. At 
the maximum size of our Monte-Carlo simulation L = 24, 
the potential is accurate to the eighth decimal place with 
only \ni\ =3 in each axis, and e = \pK . 



B. Finite-size effects 

In the vicinity of Wigner crystallization, the finite-size 
effects are very strong. The size dependence of the single 
particle density of states obtained from Monte-Carlo data 
for c? = 3, a = 0.3, and T = 0.0554 is shown in figure 4. 
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Figure 4: (Color online) (a) Two-Gaussian fit of the single 
particle density of states from different size of MC simula- 
tion. At this particular temperature, They fit perfectly with 
the EDMFT results. The fit form is shown in the inset, 
(b) The finite size scaling of the distance (left) and width- 
squared (right) of the two-Gaussian function for different sizes 
L =8,10,12,16,18,20,24. 



The nonlinear (two-Gaussian) fitting is done using 
IGOR 6.01. The fitting parameters, i.e., the distance 
between the Gaussian peaks and width squared as a func- 
tion of is shown in the Fig 4(b). This allows us to 
perform an accurate extrapolation to i = cxj, and the re- 
sult is found to be in excellent agreement with EDMFT 
prediction. Note how the finite-size result remains very 
far from the L ~ oo extrapolant even for our largest sys- 
tem size (L = 24). Accurate results, thus, simply cannot 
be obtained without such finite size scaling analysis. 



6 



1 P. W. Anderson, Phys. Rev. 109, 1492 (1958). 

2 N. F. Mott, Proc. Phys. Soc. (London) A62, 416 (1949). 
^ E. Wigner, Phys. Rev. 46, 1002 (1934). 

* B. Tanatar and D. M. Ceperley, Phys. Rev. B 39, 5005 
(1989). 

^ J. Huang, D. S. Novikov, D. C. Tsui, L. N. Pfeiffer, and 

K. W. West, Physical Review B 74, 201302 (2006). 
^ G. S. Boebinger, Y. Ando, A. Passner, T. Kimura, M. 

Okuya, J. Shimoyama, K. Kishio, K. Tamasaku, N. 

Ichikawa, and S. Uchida, Phys. Rev. Lett. 77, 5417 (1996). 

C. Panagopoulos and V. Dobrosavljevic, Phys. Rev. B 72, 

014536 (2005). 
® E. Dagotto, Science 309, 257 (2005). 

® N. F. Mott, Metal- Insulator Transition (Taylor & Francis, 

London, 1990). 
^° B. Spivak, Phys. Rev. B 64, 085317 (2001). 
" A. Kosevich, Sov. Phys. JETP 50, 1218 (1979). 
^2 A. L. Efros, Phys. Rev. Lett. 68, 2208 (1992). 

J. S. Thakur and D. Neilson, Phys. Rev. B 54, 7674 (1996) 

A. P. Ramirez, Annu. Rev. Mater. Sci. 24, 453 (1994). 

A. A. Pastor and V. Dobrosavljevic, Phys. Rev. Lett. 83, 

4642 (1999). 

" R. Chitra and G. KotUar, Phys. Rev. Lett. 84, 3678 (2000). 
S. Pankov and V. Dobrosavljevic, Phys. Rev. Lett. 94, 



046402 (2005). 

J. Schmalian, D. Pines, and B. Stojkovic, Phys. Rev. B 60, 
667 (1999). 

M. P. Lilly, J. L. Reno, J. A. Simmons, L B. Spielman, 
J. P. Eisenstein, L. N. Pfeiffer, K. W.West, E. H. Hwang, 
and S. DasSarma, Phys. Rev. Lett. 90, 056806 (2003). 
A. Georges, G. Kotliar, W. Krauth, and M. J. Rozenberg, 
Rev. Mod. Phys. 68, 13 (1996). 

P. Werner, A. Comanac, L. deMedici, M. Troyer, and A. J. 
Milhs, Phys. Rev. Lett. 97, 076405 (2006). 
P. Ewald, Ann. Phys. 369, 253 (1921). 
U. Essmann, L. Perera, M. L. Berkowitz, T. Darden, 
H. Lee, H. G. Pedersen, J. Chem.Phys. 103, 8577 (1995). 
E. R. Smith, Proc. R. Soc. London 375, 475 (1981). 
M. Miiller and S. Pankov, Phys. Rev. B 75, 144201 (2007). 
Strictly speaking, the nonlocal effects ignored by EDMFT 
are, for a <^1, negligible throughout the broad pseudogap 
regime, but not in the narrow critical regime close to Tc- 
The smallness of the melting temperature for a classical 
(continuum) Wigner crystal can similarly be understood^ 
by comparing it to the Debye temperature of sheer 
phonons. 



